Modified Becke-Johnson potential investigation of half-metallic Heusler compounds 
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We have investigated the electronic structures of various potentially half-metallic Heusler com- 
pounds with the Tran-Blaha modified Becke-Johnson (TB-mBJLDA) potential within the density 
functional theory. The half-metallic gaps are considerably enhanced with respect to values from 
the Perdew-Burke-Ernzerhof (PBE) functional. In particular the unoccupied densities of states are 
modified by mBJLDA, and agreement with experiment is considerably worse than for PBE. The 
agreement of the densities of states can be improved by reducing the Tran-Blaha parameter c. How- 
ever, ground state properties such as the hyperfine fields are more accurately described by PBE than 
by mBJLDA. Despite its success for ionic and covalent semiconductors and insulators, we conclude 
that mBJLDA is not a suitable approximation for half-metallic Heusler compounds. 



I. INTRODUCTION 

The half-metallic Heusler compoundai"— have fre- 
quently been considered as ideal electrode materials for 
spintronic devices. High tunnel and giant magnetore- 
sistance and high spin injection efficiency are expected 
from half- metals, i.e., materials that have full spin po- 
larization at the Fermi energy. Heusler compounds are 
ternary intermetallic compounds with the general chemi- 
cal formula X2Y Z , where X and Y are transition metals, 
and Z is a main group element. They form the cubic L2i 
structure, which has inversion symmetry and belongs to 
space group Fm3m. Some of the closely related inverse 
Heusler compounds with the Hg2CuTi prototype struc- 
ture (space group F43m) have only recently been discov- 
ered to exhibit half-metallic ferromagnetism as wellj^>^ 
The half-metals from both classes of compounds follow 
the Slater-Pauling rule, which relates the magnetic mo- 
ment m (given in /ie per formula unit (f.u.)) and the 
number of valence electrons Ny via m = Ny — 24.— 

Most electronic structure studies of Heusler com- 
pounds are based on the Kohn-Sham framework of den- 
sity functional theory^i^* (DFT), which is today the main 
tool to obtain the electronic structure of solids. How- 
ever, a well-known failure of this framework is the un- 
derestimation of band gaps. This is closely related to 
a missing derivative discontinuity Axe in the approxi- 
mate exchange-correlation (xc) functionals. However, the 
Kohn-Sham gap eg differs from the true gap Eg by this 
discontinuity even for the exact xc functional, which has 
to be computed and added to the Kohn-Sham gap "by 
hand" . ~~ This is in principle also true for the gap of half- 
metals,^^ 

The appropriate framework to discuss band gaps is 
the many-body perturbation theory, e.g., within the GW 
approximation.^^ Unfortunately, this approach is com- 
putationally very expensive. Tran and Blaha recently 
proposed an alternative, equally accurate and computa- 
tionally cheaper method to obtain the gap directly as 
differences of Kohn-Sham eigenvalues: they modified the 
Becke-Johnson exchange potentiali^ with an additional 



parameter c, so that it readf 



,,mBJ/ 



(r)=<^^(r) + (3c 




(1) 



where n^(r) is the spin-dependent electron density 
and icr(r) is the spin-dependent kinetic-energy density. 
v^^{r) is the Becke-Roussel potential, which models the 
Coulomb potential created by the exchange hole.— Due 
to the kinetic-energy dependent term in the mBJ poten- 
tial, it reproduces the step-structure and derivative dis- 
continuity of the effective exact exchange potential of free 
atoms. The parameter c was proposed to be determined 
self-consistently from the density by 
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with two parameters q;,/3, which have been chosen to fit 
the band gaps of a broad range of sohds. It can be related 
to the dielectric response of the systemJ^ii^ c increases 
with the gap size and has a typical range of l.l-ljj^ 
The mBJ potential has been proposed to be combined 
with LDA correlation (mBJLDA). Its particular merits 
and limits have been reviewed by Koller et alM In a 
recent paper, Koller et al. have suggested a new and more 
balanced parametrization of c, based on a larger test set 
of solids.^'' This reparametrization gives, however, rather 
similar results as the original parametrization for small 
gap materials. Making use of the kinetic-energy density, 
mBJ is formally a meta-GGA potentiali^i 

In a recent paper, Guo and Liu have used the mBJLDA 
potential to investigate the half-metallic ferromagnetism 
of zinc blende transition metal pnictides and chalco- 
genides. They found that mBJLDA enhances the half- 
metallic gaps significantly with respect to conventional 
DFT calculations.— In the present paper, we aim to in- 
vestigate if the half-metallic gap of Heusler compounds 
is enhanced with mBJLDA, and if such an enhancement 
leads to an improved description of the electronic struc- 
ture. 
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II. COMPUTATIONAL DETAILS 

All calculations in this work are based on the full- 
potential linearized augmented-plane-wave (FLAPW) 
method. The mBJLDA calculations are done with the 
ELK code<^ The mBJ exchange potential is available 
through an interface to the LiBXC libraryi^ F-centered 
21 X 21 X 21 k-point meshes are used with 286 points 
in the irreducible wedge for Heusler compounds and 506 
points for the inverse Heusler compounds. A gaussian 
smearing of 1 mHa is applied in all calculations. The 
mufhn-tin radii are 2.0bohr, and the momentum cutoff 
for the plane- wave expansion is fc^ax = 4.0bohr~^. The 
angular momentum expansion of potential and wavefunc- 
tions is taken to /max = 10. The mBJ exchange potential 
is coupled with the Perdew-Wang LDA correlation.^^ For 
comparison with results from a generalized gradient ap- 
proximation, we choose the PBE functional, which is 
the standard functional in most current studies of Heusler 
compounds. All calculations are based on experimental 
lattice constants. 



III. RESULTS 
A. Gaps and magnetic moments 

We have chosen the Heusler compounds for our study 
along the lines of Ref. 0, and we added some (inverse) 
Heusler compounds of current interest.—^— The main re- 
sults of our calculations are collected in Table U which 
displays the magnetic moments and band gaps within 
PBE and mBJLDA, and the Tran-Blaha parameters c. 
Notably, the c parameter is always in the range 1.09- 
1.20. The corresponding density of states (DOS) plots 
are shown in Fig. [TJ 

We observe that all materials that are half-metallic 
with PBE are half-metallic in mBJLDA as well, with an 
increased gap. Some cases which have a gap in PBE 
with the Fermi energy located outside the gap (marked 
by asterisks), become half- metals in mBJLDA (C02VAI, 
CozCrGa, CozFeAl, CoaFeGa, MnaVGa). CosMnAl, 
Co2MnSn, CoaFeSi, CoaFeGe and RuaMnSb have a 
larger gap in mBJLDA, but the Fermi energy remains 
outside the gap. 

MnaVAl and MnaVGa are ferrimagnetic Heusler com- 
pounds. MnaVAl is characterized as a half-metal with 
a majority gap by PBE, and mBJLDA has only weak 
influence on the size of this gap. This goes along with 
the smallest value of c among the materials studied here. 
MnaVGa has a pseudogap around the Fermi level with 
PBE. A gap is opened by mBJLDA, and the Fermi en- 
ergy is located within the gap. FcaVAl, FcaVGa, and 
FeaTiSn, three semimetals or zero gap semiconductors 
within PBE, are predicted to be semiconductors by mBJ- 
LDA. 

NiaMnSn, CuaMnAl, and CuaMnSn were included in 
this study to observe the influence of mBJLDA on the 



TABLE I. Magnetic moments and (half-metallic) band gaps 
computed with PBE and mBJLDA (marked as mBJ) at the 
given experimental lattice constantsQ Moments are given in 
fj,B /f-u-, gaps are given in eV, lattice constants are given in 
A. The parameter c is dimensionless. Asterisks mark gaps, 
which are above or below the Fermi energy. 





A^v 














c 






C02YZ Heusler compounds 






Co2TiAl 


25 


5.85 


0.74 


1.00 


0.40 


1.00 


1.11 


1.12 


Co2TiSn 


26 


6.08 


1.96 


2.00 


0.47 


2.00 


1.16 


1.17 


C02VAI 


26 


5.72 


1.95 


2.00 


0.36* 


2.00 


0.65 


1.13 


CoaZrSn 


26 


6.25 


1.81 


2.00 


0.50 


2.00 


1.50 


1.16 


CoaCrGa 


27 


5.81 


3.01 


3.04 


0.39* 


3.00 


1.06 


1.18 


Co2MnAl 


28 


5.75 


4.04 


4.03 


0.61* 


4.04 


1.29* 


1.14 


CoaMnSi 


29 


5.65 


4.97 


5.00 


0.81 


5.00 


1.42 


1.15 


CoaMnGe 


29 


5.75 


4.93 


5.00 


0.57 


5.00 


1.49 


1.19 


CoaMnSn 


29 


5.98 


5.08 


5.03 


0.39* 


5.04 


1.36* 


1.19 


Co2FeAl 


29 


5.73 


4.96 


4.99 


0.06* 


5.00 


0.75 


1.14 


Co2FeGa 


29 


5.74 


5.04 


5.02 


0.02* 


5.00 


0.80 


1.20 


CoaFeSi 


30 


5.64 


5.97 


5.47 


0.11* 


5.79 


0.82* 


1.16 


CoaFeGe 


30 


5.74 


5.90 


5.63 


0.09* 


5.98 


0.90* 


1.20 






other Heusler compounds 






MnaVAl 


22 


5.92 


1.94 


2.00 


0.28 


2.00 


0.48 


1.09 


MnaVGa 


22 


5.91 


1.88 


1.99 


0.02* 


2.00 


0.27 


1.15 


FeaVAl 


24 


5.76 


0.00 








0.31 


1.12 


FeaVGa 


24 


5.78 


0.00 








0.39 


1.17 


r c2 J- loli 
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6.09 


0.00 
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RuaMnSb 


28 


6.20 


4.40 


4.03 


0.28* 


4.06 


0.44* 


1.19 


NiaMnSn 


31 


6.05 


4.05 


4.03 




4.17 




1.19 


CuaMnAl 


32 


5.95 


3.60 


3.51 




3.50 




1.13 


CuaMnSn 


33 


6.17 


4.11 


3.86 




3.91 




1.19 






inverse Heusler compounds 






CraCoGa 


24 


5.80 


0.35 


0.08 


0.19* 


0.03 


0.66* 


1.17 


MnaCoAI 


26 


5.84 


1.95 


2.00 


0.43 


2.00 


0.68 


1.12 


MnaCoGe 


27 


5.80 


2.99 


3.00 


0.36 


3.00 


0.76 


1.18 


FeaCoSi 


29 


5.65 


4.99 


4.96 




5.00 


0.57 


1.16 



^ Lattice constants and experimental magnetic moments are 
taken from References H 110, [H IS HlHil. 



magnetic moment and the exchange splitting of ferro- 
magnetic Heusler compounds that do not have a gap at 
all. In all three cases the magnetic moments of PBE and 
mBJLDA are close; while the agreement with experiment 
is very good for NiaMnSn, both approximations underes- 
timate the moment of CuaMnAl and CuaMnSn. This is 
surprising, seeing that mBJLDA predicts too large mo- 
ments for Fe (2.49 ^b) and Ni (0.74 /xb) and too large 
exchange splittings at the same timeJ^ Notably, PBE 
(and LDA) already predict too large exchange splittings, 
although the moments are accuratcj^l The small influ- 
ence of mBJLDA on the magnetic moment in CuaMnAl 
and CuaMnSn may be associated with the localized char- 
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FIG. 1. Density of states plots for selected Heusler and inverse Heusler compounds. Shaded blue areas correspond to PBE 
calculations, solid red lines represent mBJLDA calculations. Spin majority states are positive, minority states are negative. 
The energy scale is given with respect to the Fermi energy. 



acter of the magnetic moment.— Thus, an increase of the 
exchange splitting does not lead to an enhanced magnetic 
moment. In contrast, Fe and Ni have more itinerant char- 
acter, so the moments depend strongly on the magnitude 
of the exchange sphtting. 

The inverse Heusler compound Cr2CoGa has been pre- 
dicted as a nearly fully compensated ferrimagnet<^ It has 
a small gap slightly above the Fermi energy. In the mB- 
JLDA description, this gap is further enlarged and the 
Fermi energy moves closer to the gap, thereby further 
reducing the magnetic moment towards zero. Experi- 
mentally, a non-zero moment is observed;^ which may 
be due to atomic disorder. Mn2CoAl is a ferrimagnet, 
and it is predicted to be a spin gapless semiconductor— 
by PBE. This prediction has recently been confirmed 
experimentally. ■^'^ This is also the case within mBJLDA, 
and the minority spin gap is only slightly enhanced (sim- 
ilar to the case of Mn2VAl). Mn2CoGe is a half- metal 
in both approximations, but mBJLDA considerably en- 



hances the minority gap. Fe2CoSi is described by PBE 
as a conventional ferromagnet with a pseudogap in the 
minority spin channel around the Fermi energy. mBJ- 
LDA opens a sizeable gap and predicts Fe2CoSi to be a 
half-metal. 

We note that the magnetic moments inside the muffin- 
tin spheres are increased with mBJLDA in all cases, 
which is compensated by an antiparallel interstitial mo- 
ment (for ferromagnets) or by the antiparallel alignment 
of the muffin-tin moments (for ferrimagnets) . The inte- 
ger magnetic moments are protected by the half-metallic 
gaps. 



B. Densities of states 

In Figured] we compare the DOS from PBE and mBJ- 
LDA calculations for some selected compounds. We ob- 
serve that the effect of mBJLDA (compared to PBE) is 
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FIG. 2. Top: Dependence of the half-metallic and semicon- 
ducting gaps on the mBJ model parameter c. Co2FeSi (open 
symbols) is not a half-metal for the plotted range of c. Bot- 
tom: DOS plots for Co2MnSi with the PBE and BJLDA 
(c = 1) approximations. 



very material-dependent and non-trivial. In all cases we 
observe that the energy range of the occupied d states is 
compressed; the d band minima are raised and the states 
close to the Fermi energy are lowered. The exchange- 
splittings of the occupied d states are enhanced in all 
cases. The low-lying s states from the sp-element are 
shifted up in some cases, or remain at the PBE position. 
We further note that the effects on spin majority and mi- 
nority states are quite different. The occupied minority 
states are least affected in most cases, while the occupied 
majority states and the unoccupied minority states show 
somewhat larger changes. 

The enhancement of the half- metallic gaps is visible for 
all compounds. Mn2VAl shows only little changes with 
mBJLDA compared to PBE, which is partly due to the 
low value of c. In the case of Co2FeSi we see that the 
Fermi energy is located slightly above the bottom of the 
minority conduction band minimum in both cases. 



C. Dependence on c 

To estimate the influence of the mBJ model parameter 
c, we calculated the gaps of Co2TiSn, Co2MnSi, Co2FeSi, 
Mn2VAl, Fe2VAl, and Mn2CoGe with different fixed val- 
ues of c. The results are displayed in Figure [2] (top). We 
can identify two classes of dependencies on c within the 
range of our investigation: an approximately linear one 
(Fig. [5] top left), and one that saturates at rather low 



c values (Fig. [2] top right). Remarkably, the three com- 
pounds belonging to the second class are ferrimagnets, 
whereas the compounds in the first class are ferromag- 
nets or nonmagnetic (Fe2VAl). In contrast to semicon- 
ductor and insulator gapsp^ the half-metallic gap does 
not necessarily grow monotonously as c is increased. 

Setting c « 0.95 nearly restores the PBE values for 
the gaps in most cases. Also the densities of states agree 
very well with those from PBE with c = 0.95. This agrees 
with the finding by Koller et al., that the magnetic mo- 
ment of Fe can be tuned to the PBE value with c slightly 
lower than unity. However, they also mention that the 
calculated hyperfine fields at the Fe sites are much worse 
with a small c, which indicates subtle differences in the 
description of the electronic structures, in particular of 
the s states. The original Becke- Johnson exchange po- 
tential (i.e., c = I) gives band structures very similar 
to those obtained with PBE, but with a slightly larger 
gap. A corresponding DOS plot is shown exemplarily for 
Co2MnSi in Fig. [2] (bottom). This agrees with earlier 
calculations with the BJLDA potential of semiconductor 
and insulator gaps.'^^ 

It is well-known that the local xc functionals (both 
LDA and GGA) fail to predict the ground state mag- 
netic moment of Co2FeSi correctlyi^ mBJLDA does not 
resolve this problem with the c value from Eq. [S] With 
c = 1.35 a just half-metallic ground-state can be obtained 
with the Fermi energy at the minority conduction band 
minimum, and the magnetic moment is G/ib, in accor- 
dance with experiments.'^^ As we discuss in the following, 
such a high value of c is rather unrealistic. 



D. Comparison to experiments 

The experimental determination of the half-metallic 
gap is notoriously difficult. Only indirect evidence for the 
presence and the size of such gaps is available. This evi- 
dence comes from tunnel spectroscopy of magnetic tunnel 
junctions with amorphous barriers^^— or from x-ray ab- 
sorption spectroscopy.'^^''^^ For Co2MnSi, a combination 
of the minority spin flip gap from tunnel spectroscopy of 
0.25 — 0.35 eV and the position of the x-ray absorption 
maximum of Co at its L3 edge of 0.9 ± O.I eV above 
gives an upper boundary of 1.2 ± 0.1 eV for the minor- 
ity gap. Note, however, that the Co rfDOS maximum is 
about 0.3 eV above the minority conduction band mini- 
mum in the PBE and mBJLDA calculations; the gap is 
therefore smaller than 0.9 ±0.1 eV. Thus, considering the 
available experimental data, the mBJLDA value for the 
minority gap is considerably too large and the PBE value 
is slightly too small. 

Also the spectral shapes of the Co L3 x-ray absorp- 
tion spectra are not well reproduced by mBJLDA. Con- 
sider, e.g., Co2TiSn (Fig. [3]left), which has a pronounced 
double-peak structure^ arising from a pure Co Cg and a 
Co-Ti hybrid t2g state, which are separated by 0.9 eV 
in the PBE calculation. Due to the excitation process. 
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FIG. 3. Experimental Co L3 absorption spectra of Co2TiSn 
and Co2MnSi (bottom row) and corresponding rfDOS with 
PBE and mBJLDA approximations. Data taken from Refs. 
[si and S 



a screened core-hole is formed and pulls the localized Cg 
states down by 0.3 eV, so the total separation of the peaks 
is predicted to be 1.2 eV, which is close to the experi- 
mental value of 1.3 eV (see Fig. [3]). Note that the ex- 
perimental spectra are broadened due to finite lifetime 
effects. This double-peak structure is clearly visible in 
the unoccupied PBE DOS of CosTiSn in Fig. E] In con- 
trast, the mBJLDA moves the eg peak up, such that it 
overlaps with the t2g peak, and no double-peak structure 
would be visible. A similar but less pronounced struc- 
ture is also present in the absorption spectra of Co2MnSi 
(Fig. [3] right) which is reproduced by PBE but is not 
present in the mBJLDA calculation for the same reason 
as for Co2TiSn. 

The gaps of Fe2VAl and Fe2TiSn predicted by mBJ- 
LDA should easily be detectable with optical meth- 
ods and by electrical transport. However, both com- 
pounds have been characterized as semimetals by 
experiments42r— Hence, the overestimation of the 
Heusler compound band gaps by mBJLDA is not lim- 
ited to magnetic cases, where the different spin densities 
may lead to an error in the determination of c, but is also 
found for paramagnetic materials. 

While the band gap (and actually the entire band 



structure) is a property of excited states, the magnetic 
moments and hyperfine fields are ground state proper- 
ties. The total moments were already discussed above 
and good agreement with experiment was shown for PBE 
and mBJLDA in most cases. We now turn to the hyper- 
fine fields i?hf, which we calculate from the Fermi con- 
tact term, ignoring other contributions.^"^ As mentioned 
above, the BJLDA potential (c = 1) provides densities of 
states that are very similar to those of PBE. Therefore 
we compare the hyperfine contact fields computed with 
PBE, mBJLDA, and BJLDA for five Heusler compounds. 
The hyperfine fields and the corresponding site-resolved 
magnetic moments are given in Table [Hi The PBE func- 
tional provides reasonable hyperfine contact fields with 
a mean absolute error (mean error) of 64kOe (-1-44 kOe) 
for the present test set. No clear trend towards an over- 
or underestimation of the hyperfine fields is visible. The 
PBE values are similar to those calculated by Picozzi 
et al.*'^ mBJLDA and BJLDA perform less well, with 
mean absolute errors (mean errors) of 133 kOe (-1-73 kOe) 
and 135 kOe (-1-101 kOe), respectively, where the mean 
of the experimental values is — 78kOe. Since mBJLDA 
and BJLDA produce very different results, the hyperfine 
fields depend critically on the value of c. Again, no clear 
trend towards an over- or underestimation with respect to 
experiment or with respect to the other approximations 
is obvious. Although BJLDA and PBE predict similar 
magnetic moments and DOS, there are notable differ- 
ences in the wavefunctions, particularly those of the s 
states, which provide the dominant contribution to the 
transfered hyperfine contact field. PBE performs best in 
this respect. 

To give some more insight into the origin of the hyper- 
fine fields, we break them down into core and valence con- 
tributions for Co2MnSi, where the Co and Mn 3s states 
are treated as core levels (Table IIIII) . Here one finds that 
both the core and valence contributions differ strongly 
among the three potentials. It has been shown that the 
core contribution scales directly proportional to the mag- 
netic moment within the muffin-tin spheres;^ Here, we 
find on average -130kOe//iB for PBE, -68kOe//j,B for 
mBJLDA, and -43kOe/A(B with BJLDA, with little dif- 
ference between Co and Mn. However, Novak et al. have 
shown by applying a self-interaction corrected potential 
to the core states, that the core contribution should ac- 
tually be larger in magnitude than that obtained with 
PBE.-*^ Thus, the (m)BJLDA values are clearly worse 
than the PBE values. 

The mBJLDA potential was originally designed as an 
empirical means to obtain the band gap directly as dif- 
ferences of Kohn-Sham eigenvalues.''^ Thus, the overall 
band structure is not intended to be improved over other 
semilocal approximations, such as the PBE functional. In 
other words, the mBJLDA potential is neither designed 
to be a better approximation to the quasiparticle self- 
energy, nor is it designed to be a better approximation 
to the true Kohn-Sham potential, which is refiected in 
the worse description of the hyperfine fields. As shown 
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TABLE II. Hyperfine contact fields (in kOe) and site-resolved 
magnetic moments (in ^b) of diosen Heusler compounds. See 
text for mean (absolute) errors. 





PBE 


mBJLDA 


BJLDA 


exp. 






m 


Hhf 171 


Hm 


m 










CooMnSi 








Co 


-108 


1.07 


-139 1.18 


-45 


1.12 


-145^ 


Mn 


-219 


2.85 


-154 3.01 


-51 


2.91 


-336" 


Si 


26 


-0.04 


60 -0.06 


40 


-0.05 










CooMuSu 








Co 


-139 


1.00 


-200 1.17 


-83 


1.06 


-156^ 


Mn 


-224 


3.02 


-178 3.19 


-44 


3.06 


-344_^ 


Sn 


89 


-0.03 


288 -0.04 


133 


-0.03 


105^^ 








CooTiSn 








Co 


52 


1.05 


14 1.27 


106 


1.09 


21_^ 


Ti 


-82 


-0.03 


-137 -0.19 


-85 


-0.04 




Sn 


94 


0.00 


331 0.01 


155 


0.01 


82*' 








Mn2 VAl 








Mn 


-76 


1.45 


7 1.59 


26 


1.46 


-99^ 


V 


-100 


-0.84 


-314 -1.01 


-249 


-0.85 


-64^ 


Al 


-31 


-0.02 


-37 -0.03 


-31 


-0.03 


-25^^ 








Cu2MnSn 








Cu 


-234 


0.06 


-219 0.05 


-190 


0.04 


-175^ 


Mn 


-79 


3.44 


-4 3.62 


139 


3.41 




Sn 


488 


0.00 


535 0.00 


539 


0.00 


196'' 



^ Reference I44l 
^ Reference 4a 
Reference 43 



by Kresse et al., the mBJLDA band dispersions can be 
even worse than those of PBE, for example in the case 
of optical absorption spectra of Si02.— They also point 
out, that a local potential cannot simultaneously pre- 
dict the gap and the band dispersions correctly. On the 
other hand, mBJLDA seems to be beneficial for the lo- 
cal magnetic moments of strongly correlated materials^^ 
and for the oxygen K edge energy-loss near-edge struc- 
ture (ELNES) description of NiO.^^ 

IV. CONCLUSIONS 

We have calculated magnetic moments and densities of 
states of 26 (inverse) Heusler compounds with the Tran- 
Blaha modified Becke- Johnson + LDA (mBJLDA) po- 
tential. In the half-metallic cases, the gaps are much 
wider than those obtained with the PBE functional. We 
have shown for some cases (for which sufficient experi- 
mental data are available) that mBJLDA docs not im- 
prove the description of the half-metallic band gaps or 
band structures of these compounds with respect to the 
PBE functional. The original Becke- Johnson potential -I- 
LDA (BJLDA) predicts similar band structures as PBE, 
but with a slightly larger gap. Also the magnetic mo- 
ments agree well with PBE. Such a description does seem 
to be more reasonable for materials with a metallic di- 
electric response, whereas highly correlated magnetic in- 
sulators have been shown to be overall better described 
by the mBJLDA potentialji^i^ On the other hand, the 
hyperfine fields (as an important ground state property) 
predicted by mBJLDA and BJLDA are worse than those 
from PBE. 



TABLE III. Total contact hyperfine fields, valence, and core 
contributions for Co2MnSi with PBE, mBJLDA, and BJLDA. 







PBE 




mBJLDA 


BJLDA 




total 


val. 


core 


total val. core 


total val. core 


Co 


-108 


32 


-140 


-139 -61 -78 


-45 1 -46 


Mn 


-219 


152 


-371 


-154 55 -209 


-51 79 -130 


Si 


26 


28 


-2 


6 62 -2 


40 41 -1 
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